Solvent Entropy in and Coarse-Graining of Polymer Lattice Models 
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The self- and mutual-avoiding walk used in conventional lattice models for polymeric systems 
requires that all lattice sites, polymer segments, and solvent molecules (unoccupied lattice sites) 
have the same volume. This incorrectly accounts for the solvent entropy (i.e., size ratio between 
polymer segments and solvent molecules), and also limits the coarse-graining capability of such 
models, where the invariant degree of polymerization controlling the system fluctuations is too 
small (thus exaggerating the system fluctuations) compared to that in most experiments. Here 
we show how to properly account for the solvent entropy in the recently proposed lattice models 
with multiple occupancy of lattice sites [Q. Wang, Soft Matter 5, 4564 (2009)], and present a 
quantitative coarse-graining strategy that ensures both the solvent entropy and the fluctuations in 
the original systems are properly accounted for using such lattice models. Although proposed based 
on homogeneous polymer solutions, our strategy is equally applicable to inhomogeneous systems 
such as polymer brushes immersed in a small- molecule solvent. 

PACS numbers: 05.10.Ln, 61.25.he, 82.60.Lf 
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Lattice models, where the continuum space is divided 
into a lattice, are ubiquitously used in both theoret- 
ical and simulation studies of polymers; perhaps the 
most well-known example is in the original Flory-Huggins 
theory In such conventional lattice models, each lat- 
tice site is occupied by at most one polymer segment rep- 
resenting a group of real monomers, and an unoccupied 
lattice site is often treated as a solvent molecule; in accor- 
dance with this self- and mutual-avoiding walk (SMAW), 
the nearest-neighbor interactions are used. We note that 
Monte Carlo (MC) simulations using conventional lat- 
tice models are already much faster than off-lattice sim- 
ulations by avoiding the time-consuming pair-potential 
evaluation in the latter, which leads to their wide appli- 
cations. 

There are, however, some disadvantages of using a lat- 
tice model. For example, a lattice is anisotropic and 
may cause some artifacts. While this problem cannot 
be completely eliminated, one of us recently showed how 
to quantify and thus minimize the lattice anisotropy. [2j 
More significantly, in conventional lattice models, SMAW 
requires that all lattice sites, polymer segments, and sol- 
vent molecules have the same volume. Since a poly- 
mer segment here is the coarse-grained representation 
of a group of real monomers, this incorrectly accounts 
for the solvent entropy (i.e., the size ratio between poly- 
mer segments and solvent molecules). It also limits the 
coarse-graining capability of such models. For exam- 
ple, Binder and co-workers showed that one bond in 
the bond-fluctuation model0] maps to five consecutive 
united-atom bonds of polyethylene ;[J] chains of A" ~ 10 2 
segments in this lattice model, which represents in most 
cases the upper-limit of A" in many-chain simulations 
with hard excluded- volume interactions (e.g., SMAW), 
are still too short compared to those used in experiments. 

The fluctuations in polymeric systems are mainly con- 
trolled by the invariant degree of polymerization Af = 



, where n denotes the number of chains in vol- 
ume V and R e the root-mean-square chain end-to-end 
distance. Essentially, Af quantifies how many chains a 
single chain interacts with within its volume; when A/ is 
large, the system fluctuations are small and the mean- 
field approximation neglecting fluctations/corrclations 
becomes more accurate. Taking polymer melts as an ex- 
ample, we have R e — \/N — lb, where b denotes the root- 
mean-square bond length of an ideal chain, and with hard 
excluded- volume interactions V ~ nNb 3 ; this means Af is 
on the same order of magnitude as N. Since Af is at least 
10 3 in most experiments, conventional molecular simula- 
tions with hard excluded-volume interactions (thus small 
Af- or N -values) exaggerate the fluctuations in many- 
chain systems such as concentrated polymer solutions or 
melts. 

Recently, one of us proposed fast lattice Monte Carlo 
(FLMC) simulation, where multiple occupancy of lat- 
tice sites (MOLS) is allowed with a proper Boltzmann 
weight. Q The use of soft potentials that allow particle 
overlapping (e.g., MOLS), as in its off-lattice analog [6], 
greatly speeds-up the sampling of configuration space. 
More significantly, V ~ nNb 3 (for polymer melts) no 
longer holds, and AT and N are thus decoupled; in such 
simulations one can therefore study systems at much 
larger A/ (e.g., with experimentally accessible fluctua- 
tions) while N becomes a chain discretization (coarse- 
graining) parameter that does not correspond to the ac- 
tual chain length used in experiments. This point is cru- 
cial for understanding coarse-grained models with soft 
potentials. MOLS further allows the use of Kronecker 
(5-function interactions (which are isotropic on any lat- 
tice) instead of the nearest-neighbor interactions (which 
are anisotropic except in 1D).|5| As demonstrated in 
Ref. [U, FLMC simulations (with MOLS and Kronecker 
^-function interactions) are much more efficient, in the 
study of equilibrium properties of soft matter such as 
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polymers, than both conventional lattice MC simula- 
tions (with SMAW and nearest-neighbor interactions) 
and fast off-lattice MC simulations (with pair-potential 
calculations) @ . Last but not least, when compared 
with the corresponding lattice field theories based on the 
same Hamiltonian, FLMC simulations provide a power- 
ful means for unambiguously and quantitatively revealing 
the effects of long- wavelength correlations/fluctuations in 
the system, as demonstrated in our recent workjij 043- 

In this Letter, we show how to properly account for 
the solvent entropy in lattice models with MOLS, and 
how to use such models as quantitative, coarse-grained 
description of polymeric systems, thus overcoming the 
aforementioned drawback of conventional lattice models. 

Let us consider, as the original system for coarse- 
graining, a homogeneous and incompressible mixture of 
n homopolymer chains each of N m monomers and ns sol- 
vent (S) molecules. With v and v$ denoting the volume 
of each monomer and solvent molecule, respectively, the 
average volume fraction of polymers in the mixture is 
given by (f> = nN m v/V, where V — nN m v + ngwg. In 
the coarse-grained model, we have n chains each of N 
segments and ns solvent molecules on a lattice. Each 
polymer segment or solvent molecule occupies one lat- 
tice site, and each lattice site is occupied by any com- 
bination of i > polymer segments and j > solvent 
molecules such that i + j/r = po (note that r > and 
po > here do not have to be integers). The average 
volume fraction of polymers in the coarse-grained model 
is then given by 4> — nNl 3 j poV, where the lattice volume 
V = (nN/po + ns/rpo)l 3 and / is the lattice spacing. 

To quantitatively map the original system to the 
coarse-grained model, the values of r, po, and I need 
to be determined, r represents the size ratio between 
a polymer segment and a solvent molecule, and setting 
r = N m v/Nvs gives the same <p for the original and 
the coarse-grained systems. Setting p — Nl 3 /N m v fur- 
ther gives the same V. With this mapping, it is easy 
to show that our non-bonded interaction energy H E in 
the coarse-grained model (given by Eq. © below) is 
also the same as that in the original system given by 
/3U E = X 0(l - 4>)V/v s ^, where /3 = l/k B T with k B 
being the Boltzmann constant and T the thermodynamic 
temperature, and \ is the Flory-Huggins interaction pa- 
rameter between a monomer and a solvent molecule. Fi- 
nally, to obtain the same Af, we set R e to be the same, 
i.e., 

R e ,o(N m ,X,4>) = Re,CG(N,x,<t>,r,po,b), (I) 

where the subscripts "O" and "CG" denote, respectively, 
the original and the coarse-grained systems. Note that b 
is only a function of I, which depends on the lattice type: 
for example, b = I for the simple cubic lattice (SCL) and 
V2Z for the face-centered cubic lattice; for a chosen lat- 
tice, I can therefore be solved from Eq. (JTJ. Under the 
assumption of ideal chain conformations good for poly- 
mer melts ((j) — I) or solutions in a ^-solvent (x = 0.5), 



TABLE I: Brush systems studied in some neutron-reflectivity 
experiments. 



Ref. 


Brush 


Solvent 


N m 


a 


[16] 


d-PS a 


dioctyl phthalate 


1515 


12.76 


[12 


PS 


d-toluene 


1008 


27.19 



a "d" denotes deuteration, and PS refers to polystyrene. 



TABLE II: Polymer and solvent densities. 



Polymer/Solvent 


Density (g/cm 3 ) 


Ref. 


PS 


1.047 


[18] 


d-PS 


1.12 


[19] 


dioctyl phthalate 


0.984 


[20] 


d-toluene 


0.943 


[19] 



Eq. (TTJ) gives 

I = y/(N m -l)/(N -l)a(b/0 _1 ( 2 ) 

with a being the polymer statistical segment length. 

This simple yet quantitative coarse-graining strategy 
for lattice models with MOLS is the central result of this 
Letter. It ensures that both the solvent entropy and fluc- 
tuations in the original system are properly accounted for 
in the coarse-grained model, and is equally applicable to 
inhomogeneous polymeric systems. To demonstrate this, 
we apply it to inhomogeneous systems of polymer brush 
immersed in a small-molecule solvent. 

A polymer brush refers to chains end-grafted on a sub- 
strate at high densities such that the excluded- volume in- 
teractions among adjacent chains make them stretch per- 
pendicular to the substrate. [ll| Here we consider the sim- 
plest case of uncompressed homopolymer brush on a flat 
substrate, which has been widely studied by experiments, 
molecular simulations, and various theories. [§1 IsL fill Il2j 
In particular, neutron-reflectivity (NR) studies by sev- 
eral groups have measured the polymer volume fraction 
profile in the direction perpendicular to the substrate, de- 
noted by 4>{x). It was such studies[13| in early 1990s that 
unambiguously confirmed the parabolic profile of 4>(x) 
for brushes in a good solvent as predicted by the infinite- 
stretching limit (ISL) theory of Milner et al. [3] , instead 
of the step-function-like profile assumed in the scaling 
theory of Alexander and de Gennesfla]. Table fl] lists the 
brush systems studied in some experiments, where the 
dimensionless chain-grafting density a is defined as the 
number of chains per substrate area of and Table HT1 
lists the polymer and solvent densities used in our mod- 
eling. In addition, we use a = 0.67 nm for polystyrene 
(PS) regardless of deuteration. 

In our coarse-grained lattice model, the first segment 
of all chains is grafted at x = I (i.e., in the first lattice 
layer), and an impenetrable substrate is placed at x = 0, 
which cannot be occupied by polymer segments or solvent 
molecules. There are totally N lattice layers in the x- 
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direction. The canonical-ensemble configuration integral 
of this model system is 



A' 



fc=ls=lR fc s k=l r k \ fe=l 

' ^P( r )+Ps(r)/r,p I 



(3) 



where Rk,s denotes the lattice position of the s th seg- 
ment on the fc th chain, denotes the position of the k th 
solvent molecule, the summations are over all possible 
positions of all polymer segments and solvent molecules, 
respectively, and "•" means that the products before it 
do not apply to the terms after it but the summations 
before it do. Furthermore, is the Hamiltonian of the 
k th chain due to its connectivity; for single-bond lattice 
models (e.g., SCL), f3h? = if the chain connectivity is 
maintained and oo otherwise. The Hamiltonian H E due 
to non-bonded interactions is 



m E [p,ps] 



x_ 

Po 



^/S(r)p s (r), 



(4) 



where p(r) = YTk=i EfLi ^r,R fc , s and p s (v) = Y2U S r,r h 
are the microscopic number density of polymer segments 
and solvent molecules, respectively, at lattice site r. Fi- 
nally, S denotes the Kronecker (5-function. 

Eqs. ([3j and ((4]) are similar to those in our previ- 
ous work[9(, where brushes in ID with r — 1 were 
studied. Following Ref. one can derive the cor- 
responding lattice self-consistent field (LSCF) theory, 
which is a mean-field theory neglecting the system fluc- 
tuations/correlations and is exact in the limit of po — > oo 
(i.e., n — > oo, or equivalently a — > oo). With uniform 
grafting (i.e., all lattice sites at x = I are grafted with 
the same number of chains), LSCF equations for brushes 
in a good or ^-solvent need to be solved only in the 
x-direction, with five parameters determining <fi(x): N, 
(t/pq = (f>(R e /l) 2 , r, Xi an d Ao; the last one, defined as 
the fraction of the nearest-neighbor lattice sites that are 
at the same a; of a given site, is a property of the chosen 
lattice. 

We first perform LSCF calculations to examine the 
effects of solvent entropy (i.e., r), where we use N — 
N m and Ao = 2/3 corresponding to SCL. Kent et al. 
studied d-PS brushesin a ^-solvent, dioctyl phthalatc, 
where v/vs ~ 0.252. [l6j Fig. [T]shows nearly quantitative 
agreement on 4>(x) between their NR measurements and 
our LSCF calculations using the appropriate r, without 
any parameter-fitting; similar agreement is also obtained 
for other er- values studied in Ref. [l6[ (data not shown). 
We also see that r = 1, as used in conventional lattice 
models, largely overestimates the solvent entropy in this 
case, thus leading to much flatter </>(x). 

Fig. [2] shows good agreement on <p(x) between the NR 
measurements of Karim et al. on PS brush in a good sol- 
vent, d-toluene, where v/v$ ~ 0.936, [TtJ and our LSCF 
calculations with x — 0.4. The largest discrepancy is 
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FIG. 1: (Color online) Comparisons of 4>(x) obtained from 
the NR measurements by Kent et al.(l6(] and our LSCF and 
ISL calculations with the appropriate r-value; <f>(x) obtained 
from LSCF calculations with r = 1 as in conventional lattice 
models is also shown. iV = N m , \ = 0-5, and Aq = 2/3. 



0.20 



0.15 



-e- 



0.10 



0.05 



0.00 





• Exp. 


"/Ti-^ 


LSCF ! 




ISL 


l\ 
l\ 

*\ 

1 






i\ 
\ 











25 50 

x (nm) 



75 



100 



FIG. 2: (Color online) Comparisons of <j)(x) obtained from the 
NR measurements by Karim et al.JlTj and our LSCF and ISL 
calculations with the appropriate r-value. N = 7V m , \ — 0-4, 
and Ao = 2/3. 



in the depletion zone near the substrate, which was not 
found in experiments probably due to the attraction be- 
tween PS and the silicon substrate pjj- We note that, 
while conventional lattice models with N — N m can 
closely capture the solvent entropy in this system, the 
appropriate r-value must be used in coarse-grained mod- 
els with N < N m . 

Next, we use LSCF calculations on SCL 
to examine the effects of coarse-graining (i.e., 
N). Fig. [3] shows the deviation of 4>(x) ob- 
tained with N from that with N m , A^lscf = 



^ jNl(N) 



dx\ 



AlscfG?; AO - 0lscf(x; N m )} ' / Nl(N), 

vs. N for the two cases shown in Figs. [1] and [51 Clearly, 
coarse-grained models cannot exactly reproduce the 
original system in all aspects. But we see only small 
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FIG. 3: Deviation of <j>(x) obtained from LSCF calculations 
with N from that with N m , A^lscf, for the cases shown in 
Figs. [T]and [2] 



deviations within large degrees of coarse-graining, e.g., 
A^lscf < 0.01 for N > 20, indicating the success of our 
quantitative coarse-graining strategy for lattice models 
with MOLS. 

Some discussion on the quantitative comparisons with 
NR measurements is in order here. For simplicity, we 
have used Eq. @ instead of Eq. (JXJ) in the above re- 
sults. While this is less accurate for PS in d-toluene, it 
could be offset by adjusting the x-value. The unquali- 
fied polymer-substrate interaction and the fitting proce- 
dure used to obtain 4>(x) in NR measurements also pre- 
clude unambiguous comparisons with our results using a 
simplified lattice model. 

To avoid such ambiguities and to quantitatively re- 
veal the system fluctuations/correlations, we now choose 
an athermal (x — 0) system of SMAW chains with 
N m = 400 and a = 0.02R 2 e (thus $ = 0.02) on SCL 
with an impenetrable grafting substrate as the original 
system. We have i? e ,o = (35.86 ± 0.12)Z from bulk MC 
simulations, where Iq denotes the lattice spacing of the 
original system, and solve Eq. (JXJ) for I at various N with 
Re,cc(N,r = N m /N,p = Nl 3 /N m l%,b = I) obtained 
via bulk MC simulations on SCL. The inset of Fig. fU 
shows that I monotonically increases with decreasing N. 

For N = N m (thus I — l ), Fig. [4] compares tf>(x) ob- 
tained from MC simulations (where chains are grafted as 
uniformly as possible) and ID LSCF calculations of the 
original system; the difference is due to the system fluc- 
tuations/correlations neglected by the LSCF theory (23[, 
which make the chains more stretched (i.e., 4>flmc( x ) is 
flatter than </>lscf(^)) as found in our previous work on 
brushes in an implicit, good solvent @ and ID brushes in 
an explicit solvent with r = 1@. The opposite, however, 
occurs for r = 4 and 10 on SCL, as shown also in the 
figure. 

To quantify the effects of coarse- 
graining on (j>(x), Fig. [5] shows A^flmc = 



FIG. 4: (Color online) Comparisons of <j>(x) obtained from 
FLMC simulations and ID LSCF calculations with the ap- 
propriate r- value. The error bars in ^flmc(^) with N = 100 
and 40 are smaller than the symbols and thus not shown. 
N m — 400, x — 0, and Ao = 2/3. The inset shows how I 
obtained from Eq. {TJ varies with N, where the straight line 
has a slope of —0.5876 as given by the limit of an infinitely 
long chain of self-avoiding walk [2^. 
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FIG. 5: (Color online) Deviation of tf>(x) obtained from FLMC 
simulations with N from that with N m = 400, A^flmc, for 
the case shown in Fig. [4] 



^So KN) [0flmc(x; N) - 0flmc(z; N m )f /m(N), 

vs. N; we again see only small deviations within large 
degrees of coarse-graining, e.g., A^flmc < 0.01 
for N > 50. We can further quantify the 
effects of coarse-graining on the system fluc- 
tuations/correlations, measured by A0p/c = 

^J Nl{N) dx [0lscf(z; AO - 0FLMc(a;; N)] 2 / Nl{N). 

Fig. [5] shows that, while A0p/c exhibits small variation 
with as expected for inhomogeneous systems, our 
coarse-graining strategy successfully preserves its order 
of mag nitude (10~ 3 ) for N £ [3,N m \; note that, for 
small N (< 10), A0p/c is expected to decrease with 



5 



decreasing N and vanish when N = 1. 

Finally, we note that the correct solvent entropy can 
be readily incorporated in ISL theories, which, due 
to their simplicity, have been widely used for polymer 
brushes. Such theories are based on Semenov's classi- 
cal theory]!!) that neglects fluctuations around the most 
probable chain trajectory and are strictly valid in the 
limit of infinite chain-stretching. Here we consider the 
ISL theory of Amoskov and Pryamitsyn (25[ [26j], which 
is based on the Flory-Huggins theory [J and for brushes 
of finitely extensible chains (e.g., lattice polymers) in an 
explicit solvent. [13] To incorporate the correct solvent 
entropy in this theory, the right-hand-side of Eq. (4) in 
Ref. [26| needs to be multiplied by r. The predicted 4>(x) 
of this modified theory is then shown in Figs. [1] and [2j We 
see that the agreement between the ISL and LSCF results 
becomes better for higher a and smaller \ (which make 
chains more stretched), and the differences between them 
are mainly due to the classical approximation and the 



infinite-stretching assumption used in the former (which 
are well understood in the case of continuous Gaussian 
chains[28]). We also note that in ISL theories all chains 
are grafted at x — with a reflecting substrate placed at 
x = 0, unlike in our LSCF calculations; this leads to the 
difference in <fi{x) near the substrate as shown in Figs. [1] 
and[2 

To summarize, we have showed how to properly ac- 
count for the solvent entropy in polymer lattice models 
with multiple occupancy of lattice sites @, and presented 
a quantitative coarse-graining strategy that ensures both 
the solvent entropy and the fluctuations in the original 
systems are properly accounted for using such lattice 
models. Although proposed based on homogeneous poly- 
mer solutions, our strategy is equally applicable to inho- 
mogeneous systems such as polymer brushes immersed 
in a small- molecule solvent. 

This work was supported by NSF CAREER Award 
CBET-0847016. 
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